Recurrence and Polya number of quantum walks 
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We analyze the recurrence probability (Polya number) for d-dimensional unbiased quantum walks. 
A sufficient condition for a quantum walk to be recurrent is derived. As a by-product we find a 
simple criterion for localisation of quantum walks. In contrast to classical walks, where the Polya 
number is characteristic for the given dimension, the recurrence probability of a quantum walk 
depends in general on the topology of the walk, choice of the coin and the initial state. This allows 
to change the character of the quantum walk from recurrent to transient by altering the initial state. 

PACS numbers: 03.67.-a,05.40.Fb,02.30.Mv 



Over the years random walks (RW) proved to be a 
useful tool in many areas in physics as well as in other 
branches of science p|. Quantum analogues of random 
walks have been proposed by Aharonov • The quantum 
walk (QW) found a promising application in quantum 
information for the construction of fast search algorithms 
Q, which initiated considerable effort to understand all 
aspects of QWs Recently, localisation has been found 
in 2-D [i,[101 and in 1-D for a generalized QW 

Recurrence in the dynamics of physical systems is an 
important phenomenon with many far reaching conse- 
quences ^]. In the strict sense it means the periodicity 
of the system and in a broader sense it is the recurrence 
of some particular property [lo| . Recurrence phenomena 
have been studied extensively for a variety of quantum 
systems [ll[ . In accordance with the terminology for the 
classical random walks we will consider the recurrence 
probability that the walker returns to the origin in QWs. 

In the present paper we extend the concept of Polya 
number characterizing the recurrence of RW to the quan- 
tum domain. We point out the fundamental difference 
between the recurrence behaviour of classical and quan- 
tum random walks. In particular, the recurrence of QWs 
is not solely determined by the dimensionality of the 
walk, but may depend on the topology of the walk, choice 
of the coin governing the time evolution and the initial 
coin state. Moreover, due to the intimate connection be- 
tween localisation and recurrence, we find a simple crite- 
rion for the former in QWs. 

Let us consider an unbiased random walk on an infi- 
nite d-dimensional lattice starting localized at the origin 
0. The probability that the walker returns to the origin 
during the time evolution is called the Polya number of 
the walk. If the Polya number equals one, then the walk 
is called recurrent, otherwise there is a non-zero proba- 
bility that the walker never returns to its starting point. 
Such walks are transient. Polya I4I proved that one 
and two-dimensional walks are recurrent while for higher 
dimensions the RWs are transient and a unique Polya 



number is associated to them in each dimension [13[ . 

The Polya number of a classical random walk can be 
expressed with the probability Po{t) that the walker re- 
turns to the origin in the t-ih step and its behaviour 
depends solely on the infinite sum X]t^o^'o(^)- ^ sim- 
ple criterion exists for a RW to be transient: the series 
must converge, otherwise the walk is recurrent 14]. The 
sequence Po(t) consists of positive numbers, therefore its 
decay determines the convergence of the series, if, e.g., it 
decays faster than t^^ it is convergent. 

The definition of the Polya number can be consistently 
extended for QWs by the expression 



P =1-X{{1 -Po{t)). 



(1) 



for an ensemble of identically prepared systems, in the 
following sense. Take a system and measure the posi- 
tion of the walker after one time step at the origin, then 
discard the system. Take a second, identically prepared 
system and let it evolve for two time steps, measure at 
the origin, then discard the system. Continue similarly 
for arbitrarily long evolution time. The probability that 
the walker is found at the origin in a single series of such 
measurement records is the Polya number ([T]) . The above 
definition ensures that the same criterion applies as in the 
classical case: if the series Y^'^oPoit) is convergent, then 
the walk is transient, otherwise it is recurrent Since 
any interaction with a quantum system and especially a 
measurement unavoidably disturbs it, it is important to 
fix how measurements are performed on quantum sys- 
tems when comparing to classical results (cf. definitions 



Til . l3])- One can imag- 



for hitting time in graph QWs 
ine various experimental situations where quantum walks 
are observed during propagation. The definition of re- 
currence should then be different to reflect properly the 
physics behind. 

Let us briefly review unbiased QWs on Z*^. The 
Hilbert space of the QW is given by the tensor prod- 
uct Ti. ~ Tip ® He of the position space Hp ^ £^{'£'^) 
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with the basis formed by the states |m) where m € Z'^ 
and of the c-dimensional coin space Tic- Here c is deter- 
mined by the topology of the walk, in particular by the 
number of neighbours to which the walker can jump in 
a single step. For example, the walker can jump in pos- 
itive or negative direction in a single spatial dimension 
which leads to c = 2d. Alternatively, it can perform 
such jumps in all spatial dimension simultaneously [l8| . 
so the coin space has the dimension c — 2'^. We denote 
all displacements possible in a single step of the walk by 
the set of shift vectors {e^ji = 1, . . . , c} and define the 
orthonormal basis in Tic formed by the states \ei). We 
consider unbiased walks where the shift vectors are re- 
stricted by the condition X^^i = 0- The time evolution 
operator which propagates a QW by a single step reads 

U ^ S -{Ip^C), (2) 

where Ip denotes the unit operator on the position space 
Tip, C is the coin flip operator acting on the coin state 
and the conditional step operator S is defined by 

S' = ^|m + ei)(m|® |e,)(e,|. (3) 

i 

The state of the QW after t steps 

= ^ V.(m,t)|m) ® |e,) = C/*|V(0)) (4) 

m,z 

is given by successive application of the operator ^ on 
the initial state |i/'(0))- For the above defined walk to 
be unbiased the coin flip C must be a unitary operator 
acting on TLc with all matrix elements Cy = (ei|C|ej) of 
the same absolute value 1/ \/c. Such matrices are related 
to the Hadamard matrices [l9|. 

Due to the translational invariance of the QWs in con- 
sideration, the time evolution (jl]) is greatly simplified 
with the help of the Fourier transformation 

V;(k,<)^^V(m,i)e^(-''). (5) 

m 

Here we have defined the c-component vectors 

■0(m,t) = (V'i(m,t), V2(m,i), . . . , -^^(m, t))^. (6) 

The Fourier transformation ([5]) defines an isometry be- 
tween £2(2^) and L^{W^) where K = (-7r,7r] can be 
thought of as the phase of a unit circle in and k g W^. 
The time evolution in the Fourier picture simplifies into 

Vi(k,t) = C/(k)V;(k,t- f), [7(k) = ^(k) • C, (7) 

where the time evolution operator in the Fourier picture 
U{k) is given by the coin C and the c-dimensional di- 
agonal matrix -D(k) with jth diagonal element e""'^ *', 
determined by the topology of the QW. The time evolu- 
tion in the Fourier picture ([7]) is solved by diagonalising 



the matrix U (k) which is possible due to the fact that by 
construction this matrix is unitary. Hence it has eigen- 
values of the form — exp (ia;j(k)). We denote the 
corresponding eigenvectors as Vj(\i). With this notation 
we can write the state of the walker in the Fourier picture 
at time t in the form 

^(k,i) ^^e-^C-)* (v;(k,0),z;,(k)) «,(k). (8) 
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Performing the inverse Fourier transformation we can ob- 
tain the probability amplitudes ipivn., t). As we are inter- 
ested in the recurrence of the QWs we need the proba- 
bility that the walker returns to the origin at time t 

Po{t)=p{0,t) = U[0,t)f ■ (9) 

Moreover, in analogy with the classical problem, we con- 
sider initial states which are localized at the origin, i.e. 
-0(111,0) = for all m 7^ 0. From the definition ([5|) 
follows that under such restrictions the Fourier trans- 
formation of the initial condition 7/'(k, 0) entering ([8|) is 
constant and equals ?/'(0, 0). Hence we find the exact ex- 
pression for the probability amplitude determining the 
recurrence behaviour of the QW 

mt) = j^iAti m = J (^e-^^'')*./,(k), 

/,(k) = (7/^(0, 0),t-,(k)).«,(k). (fO) 

The recurrence behaviour of a classical random walk is 
uniquely determined by its dimensionality. In contrast, 
for quantum walks, as seen from pop . one has more free- 
dom: both the initial state -0(0,0) and the coin, repre- 
sented by the phases Wj(k) and the eigenvectors Vj{k), 
can be varied. Moreover, both Wj(k) and t'j(k) are af- 
fected by the topology of the QW through the matrix 
D(k). In the following we demonstrate that the recur- 
rence of QWs can be in fact altered by exploiting the 
additional freedom offered by quantum mechanics. 

The recurrence of the QW is determined by the asymp- 
totics of (fTO)) which can be calculated e.g. by the method 
of stationary phase ji^l . Accordingly, the asymptotic be- 
haviour of Ij{t) is only affected by the points where all 
derivatives of u;j(k) vanish (saddle points). In particular, 
it is determined by the degeneracy of the saddle points, 
given by the rank of the Hessian matrix, and the cardi- 
nality of the set of saddle points. The classical result [20l | 
shows that if the phase ujj (k) has no saddle points then 
Ij{t) decays exponentially. If i^j(k) has finitely many 
saddle points which are non-degenerate, i.e. the second 
derivative of ujj (k) is non- vanishing, then asymptotically 
Ij{t) is given by Ij{t) ~ t~'^/'^. Here we assume that 
/j(k) is smooth and non- vanishing at the saddle points. 
However, if the initial state -0(0, 0) is orthogonal to some 
eigenvector Wj(k) then /j(k) = and the corresponding 
integral is zero. If this is true only for k corresponding to 
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one of the saddle points then this particular saddle point 
will not contribute to the integral Ij (t) . Hence the choice 
of the initial state might change the behaviour of the QW 
from recurrent to transient. Moreover, we can encounter 
situations where the phase Ci;j(k) does not depend ex- 
plicitly on some of the variables ki or has a continuum of 
saddle points. In such a case we can expect a slow-down 
of the decay of Ij (t) . 

Let us illustrate these features on several examples. 
We consider QWs for which the shift vectors have all 
entries equal to ±1, i.e. the coin space has the dimension 
c — 2'^. For this particular type of QWs we find that the 
diagonal matrix I?(k) can be written as a tensor product 



D{k) ^ D{ki)(g> ...(g) D{kd) 



(11) 
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of 2 X 2 diagonal matrices D{kj) — diag(e" 

First, we present examples where the phases Wj(k) 
have finitely many non-degenerate saddle points. As fol- 
lows from the above discussion the probability Po{t) be- 
haves asymptotically like t~'^ where d is the dimension of 
the walk. We start with general unbiased 1-D QW 
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-If} 



(12) 



Here the matrix U{a, /3, k) = D{k) • C(a, (3) has eigenval- 
ues Aj(a,fc) = ±e=t*"("''=) with the phase 



sina;(Q:, k) 



sin (fc — a) 

7^' 



(13) 



We find that the phase Lo{a, k) has saddle points ^ 
a ± 7r/2 and hence Po{t) behaves asymptotically like 
. Moreover, the asymptotic behaviour is indepen- 
dent of the initial state. Indeed, no non-zero initial state 
■0(0,0) exists which is orthogonal to both eigenvectors 
for fc" = a ± 7r/2. Compa ring the decay Po{t) ^ with 
the recurrence criterion 1J| we find that all unbiased 
1-D QWs are recurrent for all initial states in concord 
with the classical result. However, none of the QWs from 
the class (fT2|) exhibits localisation, since for all of them 
the probability Po(i) converges to zero. One can achieve 
localisation in 1-D by considering generalized QW s for 
which the coin has more degrees of freedom 0, [2l[ . 

Let us now turn to 2-D QWs for which the coin fac- 
torizes into a tensor product 



C7(a,/3) = C(ai,/3i)®C(a2,/?2), 



(14) 



i.e. we have an independent coin for each spatial dimen- 
sion. From the relation (jlip we find that the matrix 
U{ki, k2) also has the form of a tensor product 



U{ki,k2)^Ui{ki)®U2{k2). 



(15) 



Therefore the eigenvalues \j (k) of (fT5|) factorize into the 
eigenvalues of C/i and U2. Hence their phases i^j(k) have 



non-degenerate saddle points and we find that Po{t) be- 
haves asymptotically like . Similarly to the 1-D case 
this scaling is identical for all initial coin states. Hence 
we find that all unbiased 2-D QWs driven by indepen- 
dent coins for each spatial dimension (jl4p are transient 
for all initial states. The Polya number of this class of 
QWs is approximately 0.29143. We note that the con- 
cept of QWs with independent coins can be extended 
to arbitrary dimension d and their Polya numbers are 
uniquely determined by d [2]| . In contrast with the clas- 
sical RWs they are recurrent only for d = 1. Nevertheless, 
the transience of d-dimensional QWs with independent 
coins for d > 1 is not surprising when we recall that the 
QWs spread quadratically faster compared to the clas- 
sical RWs. However, the probability Po{t)^ which deter- 
mines the recurrence of the walk, does not need to decay 
faster for the QW compared to the classical RW. In fact, 
the constructive interference can be so strong that the 
QW will exhibit locahsation [1,13, H. 

Let us now analyze the emergence of localisation, or 
more generally the slow-down of the decay of Ij (t) . Con- 
sider the situation where ujj (k) does not depend on some 
of the variables fc^, say n of them. This opens up the 
possibility that Ij{t) in (fTO|) factorizes into the product 
of time-independent and time-dependent integrals over n 
and d — n variables. If in the reduced space oi d — n di- 
mensions a finite number of non-degenerate saddle points 
are found, then one can proceed similarly to the previous 
case and find the asymptotic behaviour 



Po{t) 



(16) 



Comparing this expression with the recurrence condition 
14| of random walks we find a sufficient condition for the 
recurrence of QWs: at least one eigenvalue of the matrix 
C(k) should be such that its phase Wj (k) depends explic- 
itly at most on a single variable and has a finite number 
of saddle points in the reduced space. Moreover, if i^j(k) 
is independent of k the leading order term of the prob- 
ability Po{t) is a constant. In the latter case Po{t) has a 
non-vanishing limit value and therefore the QW exhibits 
localisation for initial states being non-orthogonal to the 
corresponding eigenvector. However, as we have already 
pointed out, a given initial state ?/'(0,0) might lead to 
a faster decay than (|16|) if it is orthogonal to all eigen- 
vectors corresponding to such eigenvalues. Orthogonality 
must hold at least at all the related saddle points. 

The above discussed situation is nicely illustrated on 
the example of the 2-D Grover walk driven by the coin 
G with the matrix elements Gm,n = ^12l,~ ^mn- This 
QW has been extensively studied 0, 0, [l3| ■ It was first 
identified numerically 0] and later proven analytically Q 
by considering the degeneracy of the eigenvalues of the 
transfer matrix that the Grover walk exhibits localisation 
except for a particular initial state 

1 



V(o,o) = Vg = 7t(1:-i>-i>i) 



(17) 
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For the 2-D Grover walk the matrix L'^(k) has two con- 
stant eigenvalues A1.2 = ±1 and two of the form 



A3,4(k) 



cos(a;(k)) = — cos fci cos A:2- (18) 



Hence the 2-D Grover walk exhibits localisation unless 
the initial state is orthogonal to the eigenvectors corre- 
sponding to A1.2 at every point k — (fci, /c2)- The analysis 
of the eigenvectors of the matrix J7(k) reveals that such a 
vector is unique and equals (|17p . Moreover, for the par- 
ticular initial state ([T7]) the probability Po{t) decays like 
since the asymptotic behaviour is determined by the 
remaining eigenvalues A3.4(k) for which the phase uj(k) 
has only non-degenerate saddle points. Therefore the 2- 
D Grover walk is recurrent for all initial states except the 
one given in (|17p . The numerical value of Polya number 
for the particular initial state p?]) equals to that of the 
2-D QWs with independent coins [2l|. As a generaliza- 
tion of the Grover walk, one can construct for arbitrary 
dimensions a QW which is recurrent 2l|, except for a 
subspace of initial coin states. Moreover, this QW ex- 
hibits localisation in even dimensions. 

Finally, we consider the situation where the phase 
LOj{k) has a continuum of saddle points, e.g. they ahgn 
on some curve 7. The previously discussed case of Wj(k) 
which does not depend on n variables can be considered 
as a particular example of this more general situation, 
since such ojj (k) obviously has a zero derivative with re- 
spect to those n variables. The case of 2-D integrals with 
curves of stationary points are treated in [23] where it 
is shown that the continuum of stationary points slows- 
down the decay of such integrals to Similar results 

can be expected for higher dimensional saddle domains, 
however, much less is known about the stationary phase 
method for that case. 

The recurrence of the 2-D Fourier walk can be analyzed 
along these lines. It is driven by the coin F with the 
matrix elements 

F„,„ = iexp[i7r(m-l)(n-l)/2]. (19) 

Here we find that all four phases ujj (k) of the time evo- 
lution operator in the Fourier picture ?7(k) have saddle 
points at fcj = 7r/4, — 37r/4 and fcj = ±7r/2 and u;i.2(k) 
have two saddle lines 71 = (fci,0) and 72 = (fci,7r) [2l|. 
Hence Po{t) behaves asymptotically like and the 2-D 
Fourier walk is recurrent, except for the subspace 



ipFia, b) = (a, 6, a, —b) 



2\a\ 



1, 



(20) 



of the initial states which are orthogonal to the eigen- 
vectors i'i_2(k) for k G 71^2- For the initial states of the 
family ([20]) we find the asymptotic behaviour po(i) ~ 
implying that the walk is transient. This corresponds to 
the absence of the central spikes in the spatial distribu- 
tion, as it was found numerically for a special case in Q . 
The value of the Polya number for the family (PD|) can be 



varied by the complex parameters a and b of the initial 
state [ill. 

To conclude, we have extended the concept of Polya 
number to QWs in order to characterize their recurrence 
properties. Our main result is that, unlike in the classical 
case, recurrence for QWs can depend not only on the 
dimensionality of the lattice but also on the topology 
of the walk, choice of the coin operator and the initial 
state. We have formulated sufficient conditions for the 
recurrence and localisation of QWs. The present study of 
the recurrence and Polya number based on the properties 
of the matrix ?7(k) can be extended to higher dimensional 
QW s where additional interesting effects can be expected 
(2l| . The explicit dependence of the QW on the coin and 
the initial coin state opens up the possibility to design 
the value of the Polya number. 
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